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We present results of Monte Carlo simulations of random bond Potts models in two dimensions, 
for different numbers of Potts states, q. We introduce a simple scheme which yields continuous self- 
dual distributions of the interactions. As expected, we find multifractal behavior of the correlation 
functions at the critical point and obtain estimates of the exponent r\ n for several moments, n, of 
the correlation functions, including typical (n — * 0), average (n = 1) and others. In addition, for 
q = 8, we find that there is only a single correlation length exponent v describing the correlation 
length away from criticality. This is numerically very close to the pure Ising value, v = 1. 
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I. INTRODUCTION 

In the theoretical study of phase transitions in disor- 
dered systems, it is neither practical nor useful to give 
a description of a single sample with a particular real- 
ization of the disorder. Rather, a detailed theory would 
yield a probability distribution of the possible results for, 
say, the free energy, obtained by sampling over all possi- 
ble realizations of the disorder. However, one generally 
argues that important quantities, such as the free energy, 
are self averaging, that is to say independent of the par- 
ticular realization of the disorder in the thermodynamic 
limit. In this case, one is generally content to calculate 
simply the average value. This theoretical result should 
agree with an experiment value, which is done on a sin- 
gle large sample, within an error that goes to zero in the 
thermodynamic limit. 

In this paper, we will be particularly concerned with 
spin- spin correlation functions. Clearly, in a random sys- 
tem, one does not have self averaging of individual cor- 
relation functions (though the structure factor, the sum 
over all correlations, may be self averaging). Nonetheless, 
the tacit assumption is generally made that all reason- 
able ways of characterizing the distribution of correla- 
tion functions will give qualitatively similar results, and, 
in particular, will give the same values for critical ex- 
ponents. We shall call such behavior of the correlation 
functions "conventional" . 

We know, however, that there are certain cases where 
this is not true. For example, in certain models for quan- 
tum phase transitions with disorderEI, distributions of 
correlation functions are so broad that average and typ- 
ical values (characterized, say, by the median) behave in 
qualitatively different ways both at and away from the 
critical point. 

In this paper we discuss the question of whether a "con- 
ventional" description of correlation functions is correct 
for a classical model in two dimensions, the random q- 
state Potts modclB. This model has been extensij/ely 
studied by both analytical and numerical approachesoo. 
It is knownt3 that the pure model has a second order 



transition when the number of Potts states, g, is less 
than or equal to 4 and is first order for q > 4. The 
specific heat exponent of the pure system, is positive for 
q > 2, which, from the Harris criteriontj, leads to the 
conclusion that disorder will be a relevant perturbation 
in this case. Furthermorp-the. transition is always second 
order in two-dimensionalZrt3 so disorder has a particu- 
larly strong effect for q > 4, even changing the order of 
the transition. The resulting critical behavior is poorly 
understood, and its study is one of the main aims of the 
present work. 

Ludwigp has argued that the power law decay of the 
correlations at the critical point is not governed by a sin- 
gle exponent, r/, as would be expected if "conventional" 
behavior occurred, but rather shows "multifractaF'cj be- 
havior in which there is no simple relation between the 
exponents, r) n , of the various moments of the correlation 
function defined in Eq. ( |l0| ) below. Rather, the whole 
set of exponents can be conveniently represented by a 
function f(a), related to the Legendre transform of rj n , 
as we shall discuss in § [II . The r\ n have been calculated 
analyticallyu to one-loop order for q — 2 small, where the 
disorder is only weakly relevant. Away from the criti- 
cal point, Ludwigp also claims, to one-loop order in an 
expansion away from q — 2, that there is only a single ex- 
ponent v, characterizing the divergence of the correlation 
length. 

Numerical evidence that the correlations at criticality 
are multifractal has been provided by the transfer ma- 
trix calculations of Jacobsen and Cardyu. They showed 
that cumulants of the log of the correlation function stay 
non-zero as the lattice size is increased, which indicates 
multifractality. However, their values for rj n for n > 
are not obtained directly but are derived from the cumu- 
lant expansion r&pd so become inaccurate for large n. In 
addition, PicccO has argued that, for the distribution of 
disorder used in Ref. [5], the data are in a crossover region 
between random and pure behavior, at least for the larger 
q values, and so do not give accurate values of the true 
exponents of the random system. To our knowledge, the 
perturbative claim that there is only a single correlation 
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length exponent v has not been verified by other means. 

Here we study the g-state Potts ferromagnet by; Monte 
Carlo simulations, using the Wolff algorithmell which 
greatly accelerates equilibration. The main features of 
our work are as follows. 

1. We directly determine the r\ n for a wide range of 
values of q and for several values of n, including n — 
(typical), 1 (average) and others, as well as the 
whole distribution of correlation functions. We find 
that rjQ varies much more strongly with q than does 
rji- Curiously, we find that within the errorbars r]2 
varies only slightly, if at all, with respect to q. 

2. We use a different distribution for the disorder from 
that which is generally taken. It has the advan- 
tage that it is still self dualt^l, so the critical point 
can be determined exactly, but it is continuous and 
very broad. As a result, we Eucgue that it is less 
susceptible to crossover effectst3'E3 than the "bi- 
modal" distributions used before. 

3. We verify explicitly that, at least for q — 8, the ex- 
ponents for the divergence of the average and typ- 
ical correlation lengths are equal, to within fairly 
small numerical errors. Furthermore, to within er- 
rors, the value is the same as that of the pure two- 
dimensional Ising model. Note, though, that other 
aspects of the critical behavior, such as the decay of 
the correlations at the critical point, are quite dif- 
ferent from those of the pure two-dimensional Ising 
model. 

The layout of the paper is as follows. In §[n] we discuss 
the model and some characteristics of the simulations. In 
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§111 we discuss multifractal behavior, which is expected 
to describe the correlations at the critical point. Our 
results at the critical point are discussed in §IV while our 
results away from the critical point are given in §0. In 
§[v| we summarize our results and give some perspectives 
for future work. 



If K is large (low temperature) then K* is small (high 
temperature) and vice-versa. Assuming that there is a 
single transition, this must be at the self-dual point where 
K* = K = K c , with e K " = 1 + jq. 

For the random case, the model is still self dual, and 
hence at its critical point, provided that the distribution 
of the dual couplings P* is equal to the distribution of the 
original couplings. One simple example, which has been 
extensively used in numerical work, is two delta functions 
with equal weight, 



P{K) = \[6{K 



ifi) + 5(K - K 2 ) 



(3) 



which is self dual if K 2 = K\. Hence Eq. (|^) describes 
a family of self dual distributions characterized by a sin- 
gle parameter, R = K\/K%. However, if this ratio is too 
close to unity, then very large sizes are neededlij other- 
wise the system is in a "crossover" region between the 
critical behavior of the pure and the random systems. 
On the other hand, if the ratio is made too large, the dis- 
tribution is fairly close to that of the percolation problem 
at criticality, and again the system will be in a crossover 
regime. 

To avoid these crossovers, it is useful to study the 
model in Eq. (|l]) with a self-dual continuous distribu- 
tion. Although Jauslin and SwendserO describe how one 
particular continuous self-dual distribution can be con- 
structed, we are not aware of any results obtained with 
such a distribution. Here, we show that with a simple 
change of variables one can trivially generate any self- 
dual distribution. In terms of the variable y, defined by 
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the duality condition takes the simple form 
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II. THE MODEL 



The Hamiltonian of the q-state Potts model is given 



m = 



(1) 



where each site i on an N = L x L square lattice is in one 
of q-states, characterized by an integer ni — 1, 2, • • • , q. 
The couplings, K^, are positive, and include the factor 
of (3 = 1/ksT. They are independent random variables, 
drawn from a probability distribution, P(K). 

For the pure system, the partition function with cou- 
pling K is closely related to the. partition function with 
the "dual" coupling, K*, whereSB 



Hence, expressed in terms of y, any even distribution, 
Py(y), is self-dual. Note that because Eq. (|^) is so simple 
the Jacobian in the transformation from the distribution 
of y to the distribution of y* is unity, unlike the situation 
going from the distribution of K to that of K* , where 
the Jacobian is non-trivial. 

The quantity that enters in statistical mechanics av- 
erages is x = e~ K , which, for the ferromagnetic cou- 
plings discussed here, takes values in the range from 
to 1. In order that the model has strong disorder we 
seek a distribution of x that is non-zero everywhere, and 
has a finite weight at the end points. This means that 
Py(y) oc exp(— \y\) for y — > ±oo. A convenient choice, 
which we use in the rest of this paper for simulations at 
criticality, is 



Py(y) = -sech(y), 

7T 



(6) 
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for which the corresponding distribution for x = e is 



Px(x) = - 



7r (1 — x) 2 + qx 2 



(J) 



This is plotted in Fig. [l] for several values of q. Another 
advantage of this distribution is that it is easy to gener- 
ate random numbers with probability Px(x) (which are 
needed for the Wolff algorithm^) . If r is a random num- 
ber with a uniform distribution between and 1, one 
simply takes x to be 



(8) 



1 + ^/gtan(7rr/2) ' 



We apply periodic boundary conditions. The numbers 
of samples used are shown in Table |[ and the number 
of Monte Carlo sweeps are shown in Table O. Note that 
the number of sweeps for averaging is quite large, given 
that we use the Wolff algorithm which equilibrates the 
system quickly even for large sizes. The reason is that 
we wish to obtain the whole distribution of spin-spin cor- 
relation functions, so each pair correlation function has 
to be obtained with high precision. We cannot use noisy 
data for the individual correlation functions and rely on 
averaging over a large number of pairs to improve the 
statistics. This would be fine for the average but lead to 
systematic errors for other moments. For each sample we 
determine the correlation function for spins L/2 apart for 
a large number of pairs, depending on L, as indicated in 
Table EH 
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L=16 


L=32 


L=64 


L=128 


L=256 


L=512 
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100 


100 


200 


125 


18 
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3000 


1200 


400 


400 


365 


95 
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2000 


600 


210 


200 


125 


50 
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2000 


600 


200 


200 


100 
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1500 


1000 


904 


977 


426 


87 


20 


1000 


600 


513 


500 


265 


67 



TABLE I. 
q and L. 



The number of samples used for each value of 



q 


L < 256 


L=512 


2-8 
20 


50000(500) 
250000(2500) 


100000(1000) 
250000(2500) 



TABLE II. The number of Monte Carlo sweeps for aver- 
aging. In brackets is the number of sweeps for equilibration. 
Note that for q = 8, L = 512, 43 of the 87 samples were done 
with 50000(500) sweeps. 



III. MULTIFRACTAL BEHAVIOR 



L=16 


L=32 


L=64 


L=128 


L=256 


L=512 


256 


1024 


4096 


1638 


655 


10000 



TABLE III. The number of pairs used in the calculation 
of C(L/2) for different values of L. 
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FIG. 1. A plot of the self-dual distribution used in the 
simulations at the critical point, for several values of q. It is 
given by Eq. (R) of the text. 
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in which (• • ■) denotes a thermal average for a single re- 
alization of the disorder. At the critical point, the corre- 
lations decay with a power law and we find it convenient 
to define a set of exponents r\ n by 



[C{r) 



nl 1 /n 
av 



(10) 



for n = 0, 1, 2, • • • , where averages over disorder are indi- 
cated by [•••]„. The n = value in Eq. (|Tc|) , which 
gives the behavior of a "typical" correlation function, 
is obtained as the limit n — * 0, i.e. exp[ln C(r)] av . A 
"typical" correlation function can also be defined as the 
median of the distribution, with very similar results. 

In what we are calling "conventional behavior" all the 
r] n are equal. However, according to Ludwigu, there is 
no simple relation between the r\ n for the random Potts 
model in two dimensions. Instead one has multifractal 
behavior in which the probability distribution of the C (r) 
is given by 



Consider the spin-spin correlation function, C^ , for a where 
pair of sites, i and j. If the spins are separated by a 
distance r, we shall also denote this by C(r), where 



P(a) = jVexp[-/(a)mr] 

InCfr) 

a = , 

mr 



(11) 
(12) 
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and AT is the normalization. From general considerations, 
/(a) must have a minimum at some point, ao say, and the 
value of /(ao) can be absorbed into the normalization, 
so we set /(ao) = 0. 

In the thermodynamic limit, averages can be done by 
a saddle point method, with the result that for each n, 
there is a corresponding value of a given by 



and then 



/'(a) = -n (13) 

nr) n = f(a) + an. (14) 

The err or in the saddle point calculation is of ordeiH 
l/\/lnr which falls off extremely slowly with distance. 
This will be important in the analysis of the numerical 
results in the next section. 
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FIG. 2. A sketch of the function f{a) which character- 
izes the multifractal nature of the distribution of the corre- 
lation functions. To determine the exponent rj n , describing 
the power law decay of the n-th moment of the correlations 
at the critical point, according to Eq. ([H]), locate the point 
where the slope of the curve is —n and draw the tangent at 
that point. Where the tangent intersects the horizontal axis 
is 7/ n . Clearly, 770 = ao, the location of the minimum. 

Eqs. ( |l3| ) and (|lj) imply a simple graphical relation 
between /(a) and the rj n , illustrated in Fig. ||. One lo- 
cates the point on the /(a) curve with slope — n, and 
where this intersects the horizontal axis is r\ n . Clearly 
then 770 = ao, the location of the minimum. Because the 
rj n cannot be negative (otherwise the correlations would 
grow with distance), the function /(a) will diverge as a 
approaches some non-negative value. 



IV. RESULTS AT CRITIC ALITY 

We concentrate on the correlation function between 
two spins as far apart as possible in the lattice in cither 
the horizontal or vertical direction, e.g. if one site is at 
(0,0) the other is at (0,L/2) or (L/2,0). Fig. | shows 
our r esults for [C(L/2) n ]U n for several values of n at 



criticality for the case of q = 8. The slopes are equal 
to — rj n . In the "conventional" picture, they would all 
be equal. Clearly this is not the case; rather the slopes 
change with n as expected for multifractal behavior. 
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FIG. 3. Data for the q = 8 model at criticality for sizes 
between L = 16 and 512, for n — 0, 1, 2, 3 and 4. The lines are 
fits to the data and the slopes give the values for rj n indicated. 

Fig. |] summarizes our results for the exponents rj n , 
which are also presented in numerical form in Table IV. 



Notice, from Fig. |3J, that [C ra ] a v is an m creasing func- 
tion of n, as expected-pn general groundsQ since ln[C n ] a v 
is a convex functioned of n. This implies that, for fixed 
q, r) n decreases with increasing n, as seen in Fig. ^ and 



Table |V|. Clearly 770, representing the decay of typical 
correlation functions, varies very strongly with q, while 
the other r\ n vary less strongly. For small n, r\ n increases 
as q increases, whereas for large n the converse is true. 
Within errors, 1J2 is independent of q for the range stud- 
ied. 

It would be of interest to investigate the limit q — > 00, 
but we are not able to equilibrate q values very much 
larger than 20. However, a naive extrapolation of our 
data gives lim^oo 771 = 0.37 ± 0.01. 

Note that the data in Fig. [| for q = 2 show a small 
but non-zero variation of rj n with n. However, q = 2 is 
expected to have marginal behavior, in which r) n = 1/A 
(pure value) for all n but with logarithmic corrections!], 

i.e. [C(r)"]av" ~ r -1 / 4 ln(r/ro) A " , where, for example, 
Ao = —1/8 and Ai = 0. The data cannot unambiguously 
determine the form of the logarithmic factors (though it 
is consistent with them for appropriate values of A n and 
ro). The fits used to get the data in Fig. |J therefore as- 
sume a pure power law behavior, and the apparent small 
change of rj n with n for q = 2 is presumably due to ne- 
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0.746(25) 


1 


20 


0.347(10) 


2 


20 


0.231(6) 


3 


2 


0.222(6) 


4 


2 


0.210(6) 








3 


3 


0.213(3) 


4 


3 


0.195(3) 








3 


4 


0.212(5) 


4 


4 


0.190(4) 








3 


8 


0.193(3) 


4 


8 


0.163(3) 
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20 


0.176(5) 
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20 


0.144(4) 









TABLE IV. The values of rj n , defined by Eq. (|l0|) at crit- 
icality, for different values of q and n. n = 1 corresponds 
to the average correlation function and n = to the typi- 
cal correlation function. The values for n = are obtained 
by averaging results for the exponential of the average of the 
log and the median. These agreed within expected statistical 
fluctuations. The quantity in brackets is the error in the last 
decimal place. 



gleet of the logarithmic corrections. Note in Table IV, 
that for n = 1, q = 2, (whoce the logarithmic corrections 
are predicted to be absentd) we get r\\ very close to the 
pure Ising value of 1/4, as expected. 
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FIG. 4. Results for r]„, for several values of n, plotted 
against l/q. 



FIG. 5. A comparison between our (OY) estimates of 7/1 
(relevant for the average correlation function), with those of 
Jacobsen a&d Cardyu (JC), Mente Carlo results of Chatelain 
and BercheJa (CB) and PicccO^! (P). Ref. | also has data 
from Transfer Matrix techniques which agree with their Monte 
Carlo data within the errors. 

In Fig. [5] we compare our results for r\\ with those of 
other authors. An extrapolation of results pa larger val- 
ues of q is consistent with the result-pf PiccoE3 for q = 64. 
Our results also agree with Picco'sO for q — 3. For larger 
q values they lie above those of Jacobsen and CardyO. 
This is not surprising in view of Picco'st3 claim that the 
strength of disorder used by Jacobsen and Cardy is not 
large enough for q — 8 to be in the asymptotic critical 
regime. The bimodal distribution used in those studies 
is characterized by the ratio R of the two interactions. 
Ref. H used R = 2 whereas Picco argues that R ~ 10 is 
needed. More surprisingly, our results also lie somewhat 
above those of Chatelain and BercheEl who used stronger 
disorder, R > 10 for q > 8. Perhaps in that case there is 
a crossover from the percolation critical point. Our re- 
sults do, however, also agree with the result (not shown 
on the figure) of Wiseman and DomanyE^I that, for q = 4, 
r/i = 0.290 ±0.006. For q = 3, our results also agree with 
the analytical three-loop calculation of Dotsenko et al.E3, 
though the difference between the random and pure value 
is not very great in this case. 

There are few studies of any other values for rj n with 
n 7^ |L_ Conformal field theory has been employed by 
LewisEII to obtain second order expansions for r] n ; for 
q = 3 he calculates r/o = 0.314, 772 = 0.236 (which agrees 
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with a similar work by Dotsenko et al.ca), and 773 = 0.220. 
These r-j/alues agree quite well with ours. Jacobsen and 
CardyO present results from which 770 can be inferredxB 
for q = 3 and 8. However, we find that any reasonable 
extrapolation of their data to infinite system size gives 
results which are significantly below ours, especially for 
q = 8. This may be related to Picco's remarka!3 that the 
ratio of interactions, R = 2, used-in Ref. [| is too small, 
at least for q = 8. Dotsenko et alBB perform Monte Carlo 
simulations estimating 772 = 0.227C2} for both q = 3 and 
q = 4. A more detailed analysis^] of this data yields 
0.23(2) for q = 3 and 0.229(2) for q = 4. Although these 
values are slightly smaller than ours, they support the 
trend that 772 is barely changing with q. 

The exponent characterizing the relevance of weak dis- 
order is y = otp/vp, where ap and vp are the specific 
heat and correlation length exponents of the pure sys- 
tem. This vanishes, as q — > 2 + and has the value y = 2/5 
for q = 3. LudwigQ has shown that for (n — l)y <C 1 



(15) 



For g = 3 this gives 770 — ?/i = 0.05, which agrees well 
with our result of 0.046 (with an uncertainty of around 
0.006)y|See Table IV. This is in contrast to Jacobsen and 
CardyO who stated that their data did not agree with 
this relation, though they did not quote a value for 770. 

From the values of the exponents r] n i we can construct 
the multiscaling function /(a) from. ( [l3| ) and (p~4|) . This 
is shown in Fig. |^ for q = 8. 
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FIG. 6. The solid line is a reconstruction, for q = 8, of 
the multiscaling function, f(a). It is obtained as the tan- 
gent to the dashed lines, whose formulae are 71(77,1 — a), for 
n = 0, 1, 2, 3 and 4, see Eqs. @ and (0). See also Fig. |[ 
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FIG. 7. The distribution of log 10 C{L/2), for L = 512 and 
several different values of q at criticality. 

So far, we have discussed results for the moments of 
the correlation function. Now, in Fig. (?| we show results 
for the whole distribution of log 10 C(L/2), for L = 512 
and several different values of q at criticality. The vertical 
scale is chosen so that, for a log-normal distribution, the 
tails of the curves would be two straight lines symmetric 
about the peak. The data for q = 2 is not very different 
from this, but for large q, especially q = 20, the curve is 
not only much broader but also very asymmetric. Fur- 
thermore, even in the tail to the left of peak, the data for 
q — 20 is significantly curved. It would be interesting to 
understand the behavior of the distribution in the limit 
q — > co. 

Similar data is shown in Fig. ^ for q = 8 and differ- 
ent values of L at criticality. As expected, the distribu- 
tion becomes broader for increasing L. For small L the 
data for small C(L/2) is roughly straight indicating that 
this tail follows close to a log-normal distribution. How- 
ever, for L = 512, significant curvature in the tail can be 
seen. Furthermore, attempts to fit to Eq. ( pd| ) to directly 
determine the multiscaling function f(a) (and hence to 
compare with the reconstruction of this function from 
the moments in Fig. ||) were unsuccessful. Eqs. ( |l3| ) and 
(|l4| ) and Fig. [6] depend upon a saddle point approxima- 
tion which is only valid for y/h\{L/2) 3> 1, which is not 
realizable in a Monte Carlo simulation. 

As a check on our results we also determined 77! from 
the mean square magnetization at the critical point. In 
d = 2, finite size scaling predicts 



L -20/u = L - V i 



(16) 



where (3 is the order parameter exponent, v the correla- 
tion length exponent, 
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FIG. 8. The distribution of log 10 C(L/2), for q = 8 and 
several different values of L at criticality. 



i,j \n=l / 

and p n is the fraction of sites in state n. 



q 


m 


3 


0.259(4) 


8 


0.319(7) 


20 


0.345(9) 



TABLE V. The values of 771, obtained from the mean 
square magnetization, Eq. ( |l6| ) , at criticality, for different val- 
ues of q. 

For q = 8 the data is shown in Fig. || and the values 
of 771 determined this way are given in Table The re- 
sults agree with those found from the correlation function 
C(L/2), see Table IV, to within statistical errors. 



V. RESULTS AWAY FROM CRITICALITY 

Although the distribution of correlation functions ai- 
the critical point has multifractal behavior, LudwigQ 
claims that only a single exponent v describes the di- 
vergence of the correlation length as the critical point is 
approached. 




FIG. 9. Data for [m 2 ] av , defined by Eq. (0) , at the critical 
point for q = 8. From the slope, which is —rji, see Eq. ([lil), 
we get rji = 0.319 ± 0.007. 
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FIG. 10. A scaling plot of the data for the average corre- 
lation function away from the critical point for q = 8. Here 
t, defined by Eq. (tq), is the deviation from criticality. The 
solid curve is a polynomial fit. 



There are analytical predictionsH'tilHj for how v varies 
with q for q near 2. The deviation of these results from 
the value for the pure Ising model, v = 1, is small, and 
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FIG. 11. A scaling plot of the data for the typical corre- 
lation function away from the critical point for q — 8. Here 
t, defined by Eq. (tq), is the deviation from criticality. The 
solid curve is a polynomial fit. 



numerical studies have so far been unable to resolve it. 
Here we concentrate on q — 8, where disorder changes 
the order of the transition from first to second and so the 
perturbative analytical approach is not applicable. 
We define 



t=(T— T c )/T c , 



(18) 



so the Boltzmann factors, x = e , are first calcu- 
lated at the critical point from the self-dual distribution 
in Eq. (W), and are then modified by the replacement 



x 



l/(T+f) 



away from criticality. 



According to finite size scaling, the average correlation 
function away from criticality, [C(L/2,t)] av , is related to 
that at criticality by 



[C(L/2,t)] a 
[C(L/2,0)] E 



(19) 



where / is a scaling function. There is an analogous 
expression for the typical correlation function. Figs, [l^ 
and [ll] show scaling plots, assuming this form, in which 
v has been adjusted to get the best data collapse. From 
the chi-squared of the fits we estimate 



v = 1.01 ±0.02 
v = 0.99 ± 0.02 



average 
typical, 



(20) 
(21) 



showing that the values of v agree to within the errors. 

We believe that this is the first numerical calculation 
which verifies that the exponents for the average and typ- 
ical correlation length are equal. This indicates that the 



correlations have conventional, rather than multifractal, 
behavior away from the critical point. It is also interest- 
ing that v is close to (and perhaps equal to) 1, the value 
for the pure two-dimensional Ising model. 



VI. CONCLUSIONS 

In this work we have presented results for the critical 
exponents of the random g-state Potts model, for various 
moments of the correlation functions, at and away from 
the critical point. We have confirmed in greater detail 
than before that that the correlations have multifractal 
behavior at the critical point. We have also verified, for 
the first time, that there is only a single exponent v de- 
scribing the divergence of the correlation length. 

Implicit in our discussion has been the assumption of 
universality i.e. that all the exponents 7] n , and hence 
the multifractal function, /(a), do not depend on de- 
tails of the model, such as the form of the distribution of 
disorder. While universality is well established for pure 
systems, to our knowledge, it has not been convincingly 
demonstrated in situations where there are multifractal 
correlations. 

Another issue which merits further discussion is the 
form of corrections to scaling. Generally these have a 
power law form, in which the exponent is that of the 
leading irrelevant operator. However, when there are 
multifractal correlations with a continuous spectrum of 
exponents characterized by /(a), is it possible that the 
approach to the asymptotic limit is slower? This is cer- 
tainly the case for /(a) itself, where corrections fall off 
only as 1/ vrnr, much slower than a power law. Whether 
the same slow decay of corrections also applies to the ex- 
ponents rj n is not clear to us. 

We note that in spin glasses, the numerical results, 
particularly for 77, do not seem to satisfy universality^. 
Could a resolution be that the approach to the thermo- 
dynamic limit is only logarithmic, so astronomically large 
sizes are needed to see universality? 
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